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Abstract 

We have developed methods for performing qudit quantum computation in the Jaynes-Cummings 
model with the qudits residing in a finite subspace of individual harmonic oscillator modes, res- 
onantly coupled to a spin-1/2 system. The first method determines analytical control sequences 
for the one- and two-qudit gates necessary for universal quantum computation by breaking down 
the desired unitary transformations into a series of state preparations implemented with the Law- 
Eberly scheme [1]. The second method replaces some of the analytical pulse sequences with more 
rapid numerically optimized sequences. In our third approach, we directly optimize the evolution 
of the system, without making use of any analytic techniques. While limited to smaller dimensional 
qudits, the third approach finds pulse sequences which carry out the desired gates in a time which 
is much shorter than either of the other two approaches. 
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I. INTRODUCTION 



The Jaynes-Cummings model [1, 2], describing a harmonic oscillator coupled to a spin- 
1/2 system, underlies a wide variety of potential platforms for quantum computation, such 
as atoms in cavities [1], trapped ions [3-7], superconducting circuits [2, 8-16], and clouds of 
cold atoms [17]. Because of its ubiquity, understanding how to control the Jaynes-Cummings 
model is a key step in the development of a quantum computer. Synthesizing arbitrary 
states of one [1, 6, 7, 18] or more [2, 8, 9] oscillators is a widely studied first step. Growing 
interest in quantum computation, led to several proposals to use a harmonic oscillator as a 
bus between qubits [3-5]. More recent work has focused on the controlability of the system 
[19, 20]. Here, the goal is to prove that arbitrary unitary transformations may be synthesized 
with a given set of controls, without necessarily providing an explicit algorithm to perform 
the synthesis. The most advanced experimental implementation of the Jaynes-Cummings 
model is superconducting circuits, where both state synthesis in single oscillators [10, 11], 
as well as entanglement between two oscillators [9, 12] has been studied. High fidelity qubit 
readout [13] as well as the quantum von Neumann architecture have also been demonstrated 
[14]. 

In the standard approach to quantum computing, information is stored in a series of two 
level qubits and the information is manipulated by applying one and two qubit gates. In 
most schemes, single qubit gates can be done with relatively high fidelity, but two qubit gates 
often cause problems both because control over two particle interactions is less well developed 
experimentally and because they can lead to increased coupling to the environment, leading 
to decoherence. Thus, by minimizing multiparticle interactions, higher fidelity operations 
may be possible. One approach to reducing the number of multiparticle interactions is to 
make use of c?-level systems known as qudits. By combining a few of the qubits into a 
single, larger dimensional system the gates between those qubits become manipulations of 
individual qudits, which may be more reliable. A number of issues facing qudit quantum 
computation have already been considered, such as gate decompositions [21-23], simulation 
[24], and error correction [25]. In addition, qudits may offer some advantages over qubits, 
in particular non-locality without entanglement [26] , improved detection efficiencies for Bell 
tests [27] and systems to study quantum chaos [28]. A variety of experimental systems have 
been considered as qudits, including optical systems [29, 30], superconductors [2, 8, 15, 16], 
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and atomic spins [31-34]. 

In this paper we describe the synthesis of general unitary transformations on a qudit 
system defined as a two-level system and the first n — . . . N levels of a harmonic oscillator. 
Similar studies have been reported for the circuit QED system, where the interaction between 
the super conducting qubit system and the cavity field can only be switched off by detuning 
the systems with respect to each other [15, 16]. The Jaynes-Cummings model, however, also 
describes atomic systems coupled to a quantized cavity field by a Raman process with a 
classical laser field that can be both detuned and switched completely off [1] . Harmonically 
trapped ions also implement the Jaynes-Cummings model with the possibility to resonantly 
drive an internal two-state transition, and a sideband transition, which couples the internal 
state and the motional oscillator state of the system [6]. Finally, the collective occupation 
of different internal states in an ensemble of atoms can be effectively described by oscillator 
degrees of freedom, and, e.g., by the Rydberg blockade mechanism [17, 35], one of the 
populations may be effectively limited to two values and thus implement a collective two- 
level degree of freedom in the system. These systems motivate our search for effective means 
to control the Jaynes-Cummings model, using the fast, resonant processes offered by the 
laser driven atomic systems. Our analysis uses a different approach and thus supplements 
recent work by Strauch [16], which also includes use of resonant interactions. As such, 
after introducing the basic controls in Sec. II, we show that resonant Jaynes-Cummings 
interactions are sufficient to generate arbitrary transformations on a qudit system in Sec. III. 
We introduce a semi-analytic protocol to synthesize arbitrary transformations that uses a 
combination of numerical and analytic techniques to synthesize qudit transformation much 
more rapidly in Sec IV. We also show that direct numerical optimization, without making 
use of any analytic techniques, can speed up transformation synthesis even more in Sec. V. 
Finally, in Sec. VI we extend our results to multiple modes. 



II. BASIC CONTROLS 

Our goal is to use qudits consisting of the first n = 0, . . . , N levels of a harmonic oscillator. 
To control the oscillators, we couple them to a spin-1/2 system, so that in the rotating frame 
and under the rotating wave approximation, the Hamiltonian of the Jaynes-Cummings model 
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we consider is 

H =H S + 2J H S c,k, (1) 

it 

H s = - ^A(t)a z + - X (t) (cos(0(t)K + sin(0(t))a,) , (2) 
H sc , k = l -g k {t){e^a\a_ + e^o^). (3) 

The controls are the spin's detuning, A(t), drive strength and phase, 0(*)> an d cou_ 

pling strength and phase, gk(t), (3k(t). From this point forward, we will drop the explicit 
dependence on time in the controls. These controls are available in systems such as clouds 
of Rydberg atoms [17, 35], trapped ions [6, 7] and three level atoms in cavities [1]. Here, we 
focus on the generic features available in any system with these controls. 

We begin by considering the single mode case and drop the mode index k. It is convenient 
to define the finite computational subspace, h comp , consisting of the first N + 1 states of the 
oscillators and the two spin states, represented in the following by product states {\n f), \n 1 
)}n=o...N- The system's evolution can be greatly simplified in two regimes. In the first regime 
g = 0, and the Hamiltonian 

H {1) = ~\^z + \x (cos(0K + sm(4>)a y ) (4) 

only couples states within the two dimensional subspaces hn = {\n t) , \n |)}. Let Pn = 
I n ) (n| denote the projectors onto those subspaces. Using H^ l \ we can generate an arbitrary 
rotation on the Bloch sphere, 

U {1 \e,n) = e - iena ' 2 . (5) 
Next, we choose \ = 0, so the Hamiltonian 

# (2) = -\^a z + l -g (etfato-- + e~^aa + ) . (6) 

couples the states within the two dimensional subspaces, 

/i(2 ) = |{l«- 1 t),|n4.)} n^O 
" \ |0 1) n = 0. 

We also define the associated projectors, 

p(2) = 1 1 n - 1 t ) (n - 1 t I + I n I ) (n 1 | n ^ 
' I ) (0 I I n = 



(2) 

and the set of Pauli operators on the h n subspaces with n^O, 

<y X) n =| n | ) (n - 1 | I + I n - 1 f ) (n I |, (9) 
&y,n =i\ n i ) (n - 1 1 I - i\ n - 1 1 ) (n 1 1, (10) 
a 2 ,n = - | n 4 ) (n ; I + I n - 1 1 ) (n - 1 1 1- (H) 

Defining cr^o = 0, the Hamiltonian can be written 

1 1 °° 

H m = --Aa z + ^gJ2^ {cos(P)a Xin + sm(P)a y>n ) , (12) 

n=0 

and it generates the following evolution, 



UW(g,A,p,T) = e \ / . (13) 

In the appendix, we show that it is possible to use these controls to synthesize transforma- 
tions of the form 

2V+1 

Jj(2) _ ^ e -*)m(n). ff „/2p(2)^ ^ 
n=0 

Where </>(n) and m(n) are different rotation angles and torque vectors for each subspace 
hnK The transformations Eq. (5), Eq. (13) and Eq. (14) form the basic building blocks from 
which all other controls will be built. 

III. ANALYTIC SYNTHESIS OF ARBITRARY TRANSFORMATIONS 

Our unitary design scheme builds on the state preparation protocol originally developed 
by Law and Eberly [1], which we review here for completeness. Given the controls available, 
we need to synthesize a transformation, U, which maps an arbitrary state of the oscillator- 
spin system, to |0 l). To find U, we break the problem into a series of substeps, 

N 

U-uSlWui*. (15) 

k=l 

Where Uj and have the form of Eq. (5) and (13). The pulse sequences are chosen so 
that the population is removed sequentially from each harmonic oscillator level. The state 
|0) initially has population spread over all iV oscillator levels, 

TV 

\4>) = Cn ± l n + c ™t \ n t) • ( 16 ) 

n=0 



The combination U^U^ is designed to remove all the population in h$ by transferring 
the population to \N — 1 1)- First, transfers all the population in h$ to \N \), then 
transfers all the population in hff to \N — 1 1)- The result is a new state, 

N-l 

= ^2 c «+ ' n ^ + c «t \ n t) > ( i7 ) 

n=0 

whose highest populated oscillator state is iV — 1. The rest of the pulse sequences proceed 
in a similar manner with U^U^ clearing out each oscillator level one by one until all the 
population has been transferred to A final U^ +1 transfers all of the population to |0 !). 

To see how the Law-Eberly scheme for state preparation can be extended to unitary 
transformation synthesis, we note that a unitary transformation may be defined as a trans- 
formation which maps a particular basis set \(p n ,s) back to the \n, s) basis. In other words, 
U \4> n ,s) — \n, s ), where n is the harmonic oscillator level and s = {t>4-} is t ne s P m state. 

Thus, we can break a transformation up into a series of substeps, 

f/= nri f/(ns) - 

n s 

The main goal of each substep is to complete one state transformation, which maps a given 
state \4> ns ) back to the corresponding basis state \ns), so that 

\ns) = f/ (ns) . . . U m) U (oi) \(f> ns ) . (19) 

Each £/( ns ) accomplishes this goal via a Law-Eberly type sequence of substeps. However, 
there are a couple of extra constraints on which require some modifications of the 

Law-Eberly scheme. 

After the first transformation, we will have transformed j^oj.) to |0 |). All other 

\<pns) will be transformed to a new set of target states \<pnl) = \<f> ns ). We do not want 
population leaving the computational space, so must satisfy the additional constraint 
that | (f>ns ) £ h com p. 

The next pulse sequence, will transform \4>^} to |0 f), but we also need to en- 

sure that the previously prepared state, |0 l) is unchanged. Thus, we also require that 
[/(°t) |o |) = |0 4) For all other target states will transform them to = \($). 

Once again, the transformation must prevent any population from leaking out of the com- 
putational space, so \<pns) G h comp . 
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All the other £/( ns ) will have a similar form. First of all, f/( ns ) will complete the corre- 
sponding state preparation, so that 

u {ns) i0i 2 r s+1) ) = m. (20) 

Where we use s = for spin up and s = 1 for spin down. Secondly, we must keep track 
of the changes previous U^ ns ^s have made to the original target states. The transformation 
U( ns "> takes these states from |0n 2n ~ s+1 ' ) ) to some new set of target states, |0n ™ _s+2 ' ) ), which 
must remain in the computational subspace, in order to prevent population from leaking 
out of that space, 

u( nV |^-+D) = \^ s 7 s+2) ) G h comp n'>n (21a) 
U {ni) |0£r S+1) > = |0i 2 r s+2) > G h comp n'>n-n = n' and s' =| . (21b) 

Finally, after each U^ ns \ we will have a growing set of previously prepared states, 
{|0 1) , |0 f) , |1 l) , ■ ■ ■}, and subsequent U^ ns ^s must not allow any further changes to these 
states, 

C/K) \ n 's') = \n's') ri < n (22a) 

U&V \n's') = \n's') n' <n;n = n' and s' =| . (22b) 

We define hid = {|0 I) , |0 |) , |1 4) , • • •} to be those previously prepared states that must 
remain unchanged throughout the rest of the pulse sequence. To fulfill Eq. (21)-(22) and 
carry out the requisite state preparation, all we need to do is break U ns up as we did in 
Eq. (15) 



{y-rN-n Tr (ns,2) Tr (ns,l) ± 
llfc=i u k u k ~~ I / 2 3) 

TT (ns,l) T-rN-n TT (ns,2) TT (ns,l) , 
u N-n+lilk=l U k U k b —+i 

where 

u{ ns,l) = I 4l k+hn S =| (24a) 
U JV-fc+l,n-l S = 4' 

^r ,2) =^u +1 , ra , (24b) 
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and 



z=o 

A/" 



= E ^ + <^ (2) + WSL + ^U- (25b) 



(i) 



We show in the appendix how to synthesize u^jj- and u^jj- with the available controls. 

By inspection C/ J [ ns,1) and C^ ns ' 2) above satisfy Eq. (21)-(22). They satisfy Eq. (20) if M J 
and M {2) are chosen as described in the following, and if the unitary transformations Wcomp 

(1) (2) (2) (2) (2) (2) 

and Wcomp operate within the subspaces of h comp orthogonal to (Bh\ ®h 2 ©• • -®h)J®hj 
and © © © ... © h$ © respectively, and the unitary transformations W^ omp 
and W±c 0m p operate only on the subspace orthogonal to h comp . 

We begin by considering the K th step synthesizing with s =|, as depicted in Fig. (1). 
The population of n^Ti 1 U^'^U^' 1 ^ \4>ns~ s+V> ) is spread throughout the computational 
space, with the size of the circles representing the amount of population in each level. By 
unitarity, n^Ti 1 uj^'^uj^ 18 ' 1 ^ |0n 2 s n ~ s+1 ' ) ) has no overlap with the states in h^. Fig. (la) 
depicts the action of U^ ns,1 \ It leaves unchanged while M^_ K+1 transfers all the popula- 
tion of nfji 1 U^U^ lsX) |<tf s n ~ s+1) ) within h ( H K+1 to \N - K + H). Fig. (lb) depicts the 
action of U^ s ' 2 \ which once again does not effect hid, while M^_ K+1 transfers all the popu- 
lation in h {2) _ K+1 of U { ™> 1] nfji 1 U { ™' 2) U { ™> 1] |0E n " s+1) ) to \N - K t>. The net result of the 
K th step synthesizing is shown in Fig. (lc). As we see, all the population in h$_ K+1 
°f n^Ti 1 U^^U^^ |0n 2n_s+1 ' ) ) has been transferred to h$_ K . If K ^ N — n, then we are 
only concerned with transferring population, so the torque vectors of neither M^_ K+1 nor 
K+l require a z-component. In the final step of synthesizing the transformation 

M,?^ must transfer all the population in h% of nf=T _1 U^U^ to 

\n |) and set the resulting phase in \n t) to zero. In this case, we must control the phase of 
the final state, so requires a torque vector with a ^-component. 

Similarly, If s =| and K ^ N — n + 1, then M^_ K+l K+1 must transfer all the pop- 
ulation in h$_ K+1 K . Since we are only concerned 
with transferring population, the torque vectors of neither M^_ K+1 nor M^_ K+1 require a 
^-component. During the last step synthesizing the transformation Mn^ must move 
all the population in h$ of n^i™ uj c ns ' 2 ^uj : ns ' 1 ' > |0n S n ~ s+1 ' ) } to \n \) and set the resulting 
phase in \n \) to zero. In this case, we must control the phase of the final state, so 
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(b) 



3 = n N-K + l 

(c) 

FIG. 1. The effect of U^U^ on nf^ 1 U { ™' 2) U { ™> 1] |^ 2 "" 2+1) ). The size of the circles 
represents the amount of population in each level, (a) The transformation U^ s ' 1 ^ leaves un- 
changed, while transferring all the population in h^_ K+1 of nfc=i^ uf^ 18 U^'^ \4>nP 2+1 ^) to 
\N - K + 1 1), (b) The transformation U^ s,2) again leaves hid unchanged, while transferring all 
the population in h® K+1 of Uk=i U^U^ \^>t~ 2+1) ) to \N — K f). (c) The net ef- 

feet of JJ^'^JJ^ 8 ' 1 ^ on H/fc^ l/j^ 18 U^ 18 ' 1 ^ \4>nl l 2+1 - > ) is to remove all the population from the 
highest occupied oscillator level, N — K + 1, while leaving hid unaffected. 



requires a torque vector with a z- component. 



IV. SEMI- ANALYTIC SYNTHESIS OF ARBITRARY TRANSFORMATIONS 
A. Semi-analytic protocol 

The analytic construction of the previous section is sufficient to show that the coupled 
oscillator-spin system can be controlled in a finite time but we note that the number of 
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8=1 



W 

comp 



W 

comp 



N 



(a) 

FIG. 2. Illustration of vj$ for Af+ Kn<N 

operations needed to carry out a general unitary operation on the qudit space is prohibitively 
large. In this section we will explore numerical optimization in order speed up the time 
required to synthesize an arbitrary transformation. In particular, fast implementations of 
the unitary operators u^ 1 ' 2 ^ in Eq. (25) can substantially shorten the amount of time 
required to perform general transformations. For each target transformation U, we would, 



however, have to do a separate optimization to find the corresponding set of u 



(a) 



We 



therefore proceed by synthesizing first the following transformations, which are independent 
of the chosen target transformation, 



V 



(i) 



nJ\f 



y(2) 

Vat 



V^V p(l) , mcj z /2 pW , w 



(1) 



r'(l) 



V^V p W , -i-n/2 p iira z /2 pW , VXT 



,,,nvp + W ±comp 

'(1) +W '(l) 
comp \^ ' V 



o + e~ i7r / 2 e i7T(Tz / 2 P 



(i) 

N 



w 



'(1) 



A-Comp 



±comp n = N + l,n^N 
n = H + 1, n 



N, 



3=0 



;/2p(2) , p(2) 



AT+1 



+ l y'( 2 ) +iy; (2) . 

1 comp ' ±comp 



(26a) 



(26b) 



Where w'M P and iy' (a) 



Lcomp are defined analogously to Wcolip and wj"^ and their specific 
forms are not important for the rest of the construction. Thus, V^L is the identity within 



a z-rotation within hn\ and does not couple any population out of the computational 
space. For the following constructions, we will only need n > J\f + 1. All of the transforma- 
tions V^Jj- are special cases of Eq. (14), so we have already proven we can synthesize them 
with the available controls. Once we choose the dimension of our qudit, we run a single set 
of optimizations to find the required for that dimension. 

We will need two slightly different types of pulse sequences to synthesize [T^' ^ and 



» 



jj(ns,2) ^ (jgpgjjjjjjjg on whether the torque vector of My or M^' requires a z-component. 
We begin with pulse sequences which use and to perform arbitrary rotations 



(i) 



r(2) 
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around a torque vector in the xy-plane of the hn^ and Jiff subspace, respectively, while also 
satisfying Eq. (21)-(22) 

u% ^JvWv^, (27a) 

^^VS*^. (27b) 

To understand why Eq. (27) holds, we begin by noting that if \ks) G hfj then VW) \ks) e 
h ldi so V^y/OW \ks) = VW^ lifes), and it follows that 

V$y/f$V§y/jfr) \ks) = \ks) . (28) 



Similarly, V^VO® \ns) = e - in ^^ 2 Vu^ \ns) so that 

" \ns) . (29) 
Neither nor V couple states in the computational subspace to states outside of it, 



V$Vf^V$V&&> \ns) = 



so for |0) G h comp , 



Vffiy/W) V^VO® |0) G h comp . (30) 

Finally, with M„ 2 ^ = U^ 2 \ we obtain Eq. (27a), and a similar argument yields Eq. (27b). 
The pulse sequence which synthesizes torque vectors with ^-components within h$ and 

(2) 

h n is a modification of the pulse sequence just discussed, 

u% 4^^^, (31a) 



u% =U^V^Vm V ^y/W)U^. (31b) 



If we choose 

_ e -ia(cos(v)a x +sin(u)a y )/2 ^22) 

then Mn^ = fjWljWuW has the required z-component. Similar choices for and U^' 
will give a ^-component to Mn^ = U^'^U^U^' . Finally, we note that since 
t/W" is of the form Eq. (5), we use the following slightly shorter sequence, 

«i 1 j,= y (1, V„ l SlvWV«&W". (34) 

Thus, Eq. (27), Eq. (31b) and Eq. (34) are sufficient to synthesize arbitrary unitary trans- 
formations. 
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B. Numerical optimization 



Our goal in this section is to find a pulse sequence, V^j- sa , which approximates V^j 
using fewer pulses than the analytic sequence described in Sec. III. To do so, we numerically 
optimize pulse sequences of the form 

M 

Vllsa = II Am, 0, /ft™, T ^ 2 )- ( 35 ) 



rn=l 



Where T g = 2n/g max and 

9max is the maximum allowed coupling strength. The optimization 
finds a sequence of {g^ls m -> P^l/m} which minimizes the error up to a global phase change 

1 



4(iV + l) 



mm 



Pc (V^ - e^V% a ) P C 



(36) 



V^jj- is only defined on a subspace of the computational space, and we can hence ignore 
the unimportant subspace. We begin by noting that because ^{PcV^] sa PcV^} tS aPc} < 
2(N + 1), we have 



M < 1 " yj^\^{PcV$PcV§^P c }\. 



(37) 



If we define P^X to be the projector onto the optimized subspace, then 



( N+l 

^( 2 ) + p( 2)+p ( 2+ ) i+ p 

3=1 
N+l 

Ep(2) p (2) p (2) 
r j N r N+l 

3=1 



(2) 
N+l 



N 



p(«) _ , 



AT 

Ep(2) p (2) p (2) 
3 ' N ^N+l 



(2) 
N+l 



i=i 
N 



,(2) 
JV+1 



£p(2) +p (2) +p; 

and we arrive at the following error to be minimized 



a = 1, n ^ Af+ 1, n ^ N 
a = l,n Af +1, n = N 
a = l,n = Jsf + 1, n ^ N 
a = l,n = Af+ 1, n = N 

a = 2, 



(38a) 
(38b) 
(38c) 
(38d) 
(38e) 



Ja) _ -i |.j(g) i TrfpWl/WtpWi/W p( a ) I 

_ 2(iV + 1) ±nAr iI L- r r^ l/ ? W ^nN V nN,sa^nN\- 



(39) 



Where c?^^ is the size of the subspace of the computational space orthogonal to P^. Note 
that the sums in Eq. (38) begin with j — 1 while the sums in Eq. (26) begin with j — 0. 



(a) 
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The reason is that the form of ensures that any sequence of the form Eq. (35) will leave 
|0 I) unchanged. Thus, we do not need to explicitly account for |0 I) when we perform the 
optimizations. 

For each V^fy sa) we begin with M = 3 pulses and perform at most 40 optimizations using 
MATLAB's constrained optimization routine, fmincon, with the constraint that g$m < 
g ma x- Each optimization finds a choice of {gnln, /3n,m} that is a local minimum ofe$. We 
repeat this entire procedure, incrementing M by one each time until we find a pulse sequence 
which approximates V^j- to an error better than a chosen threshold, e t hreshoid- Since our end 
goal is to find a sequence of pulses which approximates our target transformation to a desired 
accuracy, rj, and we need to apply g sa (N) = 4N 2 + 6iV + 2 transformations of the form V^j- 
to use the construction of Sec. IV A, our threshold is 

Tj 

^threshold n / , r \ • V^J 

9sa( N ) 

Finally, to bound the total time to synthesize an arbitrary transformation, we assume g max = 
Xmax, so we can bound the time for by T{U^) < \T g and T (y/ffc?) < \T g in Eq. (27) 
and Eq. (31). 

In Fig. (3) we plot the controls found by numerical optimization to synthesize V^\_i 
for N — 8. As is typical for numerically optimized controls, there is not much structure 
in gjvAr-i m- However, examination of (3^ N _ X m reveals that the torque vector lies near the 
x-axis for much of the synthesized transformation. 

Results from the optimizations are presented in Fig. (4) and Fig. (5). In Fig. (4) we 
plot the time required to reach a threshold of i] < 10~ 4 for two families of transformations, 
Vnn-i an d Vnn-2 use d to construct more arbitrary transformations. We are able to reach 
high fidelity transformations in a time less than 12T g in all cases. In Fig. (5a) we plot the 
time T{U sa ) required by the semi-analytic construction in units of T g for a target error of 
r] < 1CT 4 . Although a large number of pulses is required for the semi-analytic construction, 
it is many orders of magnitude faster than the fully analytical construction. 

Since we are able to synthesize arbitrary unitaries in dimensions up to d = 2(N + 1), we 
can compare our results to the time required to synthesize the same transformations using n 
qubits. For concreteness, we consider synthesizing an arbitrary transformation in a d — 16 
dimensional space consisting of either the N = 7 case of our oscillator and spin system or 
four conventional qubits. 
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(b) 



FIG. 3. Controls for V^ N _ 1 for N = 8. Fig. (3a) depicts the coupling strength per pulse, g$ N _ 1 m , 
while Fig. (3b) depicts the angle of the torque vector relative to the x-axis of the Bloch sphere, 



(i) 

N,N-l,m- 



The optimization finds a sequence of {g^Kj- m , m } as in Eq. (35) that minimizes 



(a) 



the resulting error, 



V. 



(a) 



V 



(a) 



nAf,sa 



The optimization is repeated a number of times with initial 



random guesses for {g^j- m , /3£fc m } and increasingly longer pulses, until a threshold error r] < 10 4 , 



|(a) 



as defined in Eq. (40), is achieved. 



The best known decomposition of an arbitrary n-qubit transformations is given in [36], 
and for 4 qubits, it requires a total of one hundred CNOTs. However, this result assumes 
any two qubits may be directly coupled while often two qubits can only be coupled via a 
bus qubit or oscillator degree of freedom. Thus for any CN0T, two extra SWAP operations are 
needed so the total number of 2-qubit operations is closer to three hundred. The number of 
two-qubit operations is a lower bound on the total time, and if each 2-qubit gate requires 
a time on the order of T g , we estimate the total time to implement an arbitrary 4-qubit 
gate with the best known gate decomposition to be bounded by T gate > 300T S . Taking this 
estimate into account, we see from Fig. (5a) that when N = 7, the semi-analytic approach 
of this section is only longer by a small factor than the bound on the multiqubit approach. 
A detailed comparison of the multiqubit and semi-analytic approach including all SWAPs, 
single qubit gates, the optimum gate decomposition in our architecture and the role of 
decoherence is beyond the present scope of this article. However, we can already conclude 
that the performance of our semi-analytic approach makes qudit computing comparable with 
qubit implementations. 

14 




FIG. 4. Time, in units of T g , required to reach a threshold of r\ < 10 4 for two families of 
transformations, V^ N _ 1 and Vjy]y_ 2 , where N is the highest controlled oscillator level 

V. FULLY NUMERICAL SYNTHESIS OF CINC' 

While the preceding section shows that arbitrary unitaries can be synthesized more 
rapidly with the semi-analytic protocol, the time required may not be optimal. To find 
a pulse sequences that performs even better, we now investigate a fully numerical optimiza- 
tion where the coupling between the spin and the oscillator is held constant while we allow 
the spin detuning, A(t), Rabi frequency, x(t), and phase, 4>(t), to be stepwise constant dur- 
ing intervals of duration dt, for a total time, Tf. The values of A(i), and <j)(t) during 
each interval are optimized with MATLAB's fminunc, which performs an unconstrained 
optimization of A(t) and <p(t). In principle, this could yield values of the Rabi fre- 
quency and detuning which are not experimentally realistic. However, we seed the initial 
random guess with x(t),A(t) < 0.9g and we find optimal controls with x(t),A(t) < 2g. 
Unconstrained optimization does not result in x{^)i A(t) 3> g because that would effectively 
decouple the spin and the oscillator, making it impossible to synthesize general unitary trans- 
formations. For each dt and Tf, 20 optimizations are performed, and the pulse sequence 
with the highest fidelity is chosen. As a test case, the gate we have optimized is 

N 

CINC' = I osc ® | | } (i | + J2 i n ® 1 ) ( n \ ® I t ) (t I- (41) 

n=0 
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FIG. 5. Time, in units of T g , to synthesize an arbitrary target transformation, U, using the semi- 
analytic construction versus the highest controlled oscillator level, N. The semi-analytic pulse 
sequence replaces the most time consuming pulses of the analytic construction with the numerically 
optimized pulses V^j- . The time plotted is the minimum time found to approximate U to the 
target error of r] < 10 -4 . 

This gate increases the oscillator level by 1, modulo N+l, conditional on the spin being 
s =j. As we will see in Sec. VI, this gate can be used to construct two qudit gates. The 
fidelity of the optimized gate, Uf n , with the target gate, U, is 

^WnTW^"'^ (42) 

For small errors, the gate fidelity can be related to the error defined above by 

F = 1 - 2r) (43) 

The optimizations have to be performed within a finite subspace and to ensure we can 
perform the optimizations in a reasonable amount of time, we keep this subspace as small 
as possible, and we penalize the leakage of population out of the computational space. In 
practice, we divide the simulated subspace into three regions given by iV < N pa d < N opt . 
The highest oscillator level in the computational space is N and the largest oscillator level 
we use in the optimization is N opt . We pad the calculation with a subspace lying above N 
but below N pad +1 and any population that leaks into this subspace is not penalized. On the 
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other hand, we penalize any population that leaks into the levels lying above N pa d during 
the course of the evolution. We can calculate the population which leaks out of N pac i during 
the course of the evolution and average over all states in the computational space, 

L = Ys dt j W W P L U fn(tj) I'- (44) 



J 

:th 



Where P is the total number of pulses, Uf n (tj) is the total evolution after the j subpulse 
and Pl is the projector onto the subspace which we wish to penalize, 

N opt 
N pad +1 

The integral over \ip) can be simplified [37], so the leakage is 

L = J2 2(N + l)(2N + 3) i^ M ^)Mm + |Tr{M(t,)}| 2 ) . (46) 

Where 

M(tj) = PcU^PlU^Pc. (47) 
Then the optimized quantity is 

C FN = l-F + wL. (48) 

After finding a pulse sequence which minimizes Cfn we calculate the system's evolution in 
a larger subspace, N c h ec k > N opt and recalculate the gate fidelity to ensure that using only a 
finite number of oscillator levels does not effect the calculation. We have found empirically 
that w = 100 ensures that the leakage out of the computational space is sufficiently small 
that when we recalculate the fidelity in the lager space, it is still high enough. 

In Fig. (6) we plot the controls for the N = 2, dt = 0.5T 9 , and Tf = 20T g case. Although 
we perform an unconstrained optimization, we see in Fig. (6a) and (6b) that neither \ nor 
A require values that are large compared to g, as anticipated. Beyond that, there is no real 
structure to the controls, as is typical when they have been numerically optimized. 

In Fig. (7a) we plot the fidelity, F versus time for iV = 2,3, with time in units of the 
vacuum Rabi period, T g . We perform the optimizations for two different sub-pulse lengths, 
dt = 0.5T g and dt = T g . We use N pad = N + 3, N opt = N + 5 and N check = 4N opt . With 
dt = 0.5T g we achieve errors on the order of 10 -4 with tens of pulses for both N = 2 and 
N = 3, substantially shorter than the semi-analytic approach, cf. Fig. (5). 
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FIG. 6. Numerically optimized control pulses used to synthesize CINC' for the N = 2, dt = 0.5T g , 
and Tf = 20T g case. 



VI. TWO MODE CONTROL 



In this section, we extend our results to cover arbitrary quantum gates between qudits 
stored in the modes of different harmonic oscillators. We consider the system described by 
Eq. (1), so that each oscillator is coupled to a central spin. In the this case, it is more 
convenient to take the computational space to consist of the first n — . . . N levels of each 
oscillator. Then if there are m oscillators, the size of the quantum comptuer's Hilbert space 
is (N + l) m . Our controls are once again the time dependent Xi 4>i A, and For 
simplicity, we assume only a single g k ^ at a time, so that only one oscillator is coupled 
to the central spin, and we optimize Xi an d A. Although the oscillators are never directly 
coupled, we can use the central spin as a bus to synthesize multiqudit gates, using the gates 
developed in previous sections. 

In analogy to qubits, multiqudit gates can be synthesized from a universal gate set. Along 
with arbitrary single qudit gates, the two-qudit CINC gate is sufficient to synthesize arbitrary 
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multiqudit gates [21]. Where CINC is defined as 

N-1 n 
CINC =^2\ ni ) (ni| ®I 2 + \N) (N\ <g) ^ | n 2 © 1 ) (n 2 |. (49) 

ni=0 ri2=0 

The action of CINC is similar to that of CINC', with the target qudit's level increased by 
one modulo N if the control is in \N). Because we have already shown we can synthesize 
arbitrary single qudit gates, showing that we can synthesize the two-qudit CINC is sufficient 
to prove we can synthesize arbitrary multiqudit gates. 

In fact, we can synthesize CINC by using the spin as a bus between the two oscillators. 
We first define, 

N-1 

BUS=J2 P n ) + ^,N + W. (50) 

n=0 

(2) (2) 

Where W acts on the h y N ' +l © h K N ' +2 © . . . subspace. Thus, assuming the central spin is 
initialized to spin down, the two-qudit CINC can be synthesized from 

BUSj 2 CINC' ls BUS s2 = CINC12 © I I ) U I + Wl- <g> I t > (t I- (51) 

Where the form of W-\ is unimportant. In the above notation M a b is meant to indicate that 
the target is system a and the control is b. Where {a, b} = {1, 2, s} is the first qudit, second 
qudit, and central spin. Thus, BUS S 2 couples the spin and the second oscillator, while CINC' ls 
couples the spin and the first oscillator. To ensure that only one oscillator couples to the 
spin we assume that gi = while we synthesize BUS S 2 and g 2 = while we synthesize CINC' ls 
We have seen in the previous section how to synthesize CINC' ls , which acts between the 
first qudit and the spin. Since BUS S 2 is a particular we can synthesize it with the 

following sequence 

M 

BUS s2 w Y[ U (2 \g max , A m , 0, dt). (52) 

m=l 

Where is taken to operate between the spin and second oscillator. As in Sec. V, we 
choose to work with a constant pulse width dt and numerically optimize A m using the same 
procedure used to synthesize V^j- sa in Sec. IV B. For each N and dt, after the optimization 
we choose the pulse sequence which achieves an error below e*- 2 -* < 10~ 4 with the fewest 
number of pulses. Finally, we synthesize BUS^ 2 with 

1 

BUS f s2 « Yl U^{g max ,-A m ,7i,dt). (53) 

m=M 
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Where the order of the product has been reversed, so that the M pulse is applied first, and 
we use the same detunings, A m used to synthesize BUS. Since A m is one of our controls, we 
are free to reverse it's sign, and (3 = n effectively reverses the sign of the coupling between 
the oscillator and the spin. 

We plot the gate fidelity vs. the total time to synthesize the two-qudit CINC in units of 
T g in Fig. (7b) for N = 2, 3 and for dt = 0.5T g ,T g . Each data point represents a different 
number of pulses used to synthesize CINC' ls , while for each N and dt, BUS S 2 is uniquely 
specified as described above. Because we have chosen to synthesize BUS S 2 with a relatively 
high fidelity, the fidelity of CINC and CINC' ls are quite close. On the other hand, BUS S 2 takes 
a finite time to synthesize, so the time to synthesize the two-qudit CINC is somewhat longer 
than the time to synthesize CINC' ls . We see that once again for both N = 2, 3 we are able to 
achieve fidelities on the order of 10 -4 . In conclusion, because the two-qudit CINC is universal 
for qudit quantum computation, we have shown that we can synthesize arbitrary multiqudit 
gates. 

VII. SUMMARY 

We have studied qudit quantum computation in the Jaynes-Cummings model. Our qudits 
consist of the first N levels of a harmonic oscillator and our controls consist of the coupling 
between the oscillator and a spin- 1/2 system as well as the drive of the spin. Beginning with 
these simple controls, we showed that arbitrary rotations within distinct two level subspaces 
are possible. Using these arbitrary rotations, along with the state preparation scheme of 
Law and Eberly [1], we designed a protocol to synthesize arbitrary single qudit gates using 
only resonant interactions between the oscillator and spin. While the analytic protocol is 
sufficient to provide a proof of principle that resonant interactions can be used to synthesize 
arbitrary transformations, it was not necessarily optimal in the time it required. However, 
we were able to use numerical optimization to reduce the time required to synthesize ar- 
bitrary transformations. A comparison between the semi-analytic routine and qubit-based 
approaches showed comparable times between the two, while our semi-analytic protocol 
had the advantage of requiring fewer multimode interactions. Furthermore, because the 
optimized interactions are relatively simple, the semi-analytic approach can be used with 
relatively large N. Although the semi-analytic approach provides a substantial speed-up 
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FIG. 7. Fidelity vs. time in units of T g for N = 2,3. In Fig. (7a), we plot the fidelity of 
CINC', which couples a single oscillator to the spin. The coupling strength between the spin and 
the oscillator is equal to the maximum Rabi frequency for the entire sequence. The spin's Rabi 
frequency, phase and detuning are piecewise constant for subpulse durations of either dt = 0.5T g or 
dt = T g and the values during each subpulse are numerically optimized. In Fig. (7b), we plot the 
fidelity of CINC, which couples two qudits stored in separate oscillators. The gate is decomposed 
into CINC' and BUS as in Eq. (51). The spin couples to one oscillator at a time, and is used as a bus 
between the two oscillators. We use CINC' as optimized in Sec. V and BUS is optimized using the 
methods of Sec. IV B. The fidelities of the two gates are so similar to one another because BUS is 
synthesized with a very high fidelity. However, a finite time is required to synthesize BUS, so CINC 
requires more time than CINC. 



over the analytic protocol, it is still not necessarily optimal. To speed up the synthesis of 
qudit transformations even further, we numerically optimized all aspects of the controls, 
without making any analytic simplifications. Because of the need to control leakage, this 
approach is limited to smaller N, but can reach high fidelities in a time an order of magni- 
tude shorter than the semi-analytic approach. Finally, we have shown that extending these 
protocols to synthesis two-mode gates with high fidelity in a short time is possible. 
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Appendix A: Synthesis of 

Our goal is to show that transformations of the form Eq. (13) are sufficient to synthesize 
more general transformations of the form Eq. (14). We cannot synthesize exactly, so 

(2) 

instead we will synthesize an approximation, Ua with an error 

E'^ = \\P C {U^ -U^)P C \\. (Al) 

Where ||M|| = A/Tr{M"l"M}. We will also bound the time needed to synthesize in 
terms of the qudit dimension, N + 1, and the desired error, E ,( - 2 \ Our method is largely 
based on [38], so we present their results for completeness. 



1. Pulse sequence 

We begin by breaking the transformation in Eq. (14) up into it's constituent Euler angles, 

= V x V 2 Vz. (A2) 

Where 

N+l 

Vl = II e iania *- n/2 , (A3a) 



?1=0 

N+l 



V 2 = Y[ e ian2ay > n/2 , (A3b) 

71=0 

N+l 

V 3 = Y[ e ia " 3a ^ /2 , (A3c) 

n=0 

are the three Euler rotations with n dependent rotation angles. By construction, = 0, 
so we are only dealing with + 1 distinct subspaces. We can perform a discrete cosine 
transform on the rotation angles [39], 

^ = X> C0S (^ N + l J' (A4) 
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The Euler rotations may then be written 



N 



V k = l[W kl . (A5) 



1=0 



Where 



N+l 



ni— I ir(n + l/2)i \ /„ 

n=0 
N+l 

nt— I ir(n + l/2)i \ ,„ 
e »a 2iv ^cos( \ + [ )a y>n /2^ 

n=0 

jV+1 

W 3l = Y\ e' Q3 ' v/ " cos ( 7r( "^+i 2) ' CTa: -")' j3: -" /2 . (A8) 



Then corresponds to the I th term in the Fourier series of the k th Euler rotation. Thus we 
have reduced the problem of synthesizing the transformations of Eq. (14), whose dependence 
on the harmonic oscillator level n is arbitrary, to the problem of synthesizing rotations around 
the x or y axis for which the n dependence is of the form y/ncos (^j^p^- 

To show that the Fourier terms may be synthesized, we begin with a sequence of pulses 
of the form Eq. (13) 

N+l N+l N+l N+l 

j- _ e ~id<f>'Vk(j x , k /2 e -id<t>'Vl<Jy,i/2 ^id<t>' y/rno x >m / 2 ^idft ^/na y>np /2 

k=0 1=0 m=0 n=0 

N+l 

« Yl e- ind(f " 7 *-" /2 . (A9) 

n=0 

Where d(f> = d(f>' 2 is small enough the above approximation holds. Repeating this procedure 
Q times allows us to generate larger rotations around the z-axis, 

N+l 

T a = t Q « Y[ e-in***,*/*. (A10) 

n=0 

Where <j>/Q — d(j) and <fi is not necessarily small. For Eq. (A10) to hold, Q must be sufficiently 
large and we will determine how large below. If we choose = 2nl/N, then we can synthesize. 
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-rr\f T ( n /7r ^ ^ t> 

W a ,kll — 1 a LJ 2 9max, U, , I ± a 

\ N + l g max J 

N+l 

e -tfij6> fciv ^cos( '■ JV+ ; )o- x , n -sin( '■ jv+ ( )a y ,„)/2^ (Alia) 

n=0 

-rr t~t f n ln d9kl ^ 
\ N + l g max J 

N+l 

~ J~J e -»rf6> fc ; v ^(cos( ^ ( ^{ 2) ' )a J; , n +sin( " (n /v l ^( 2) ' )cr yjra )/2^ (Allb) 



n=Q 

For sufficiently small, we then have 

JV+l 



w a , kll w a , kl2 » n e^"^'*^)*-.-/*. (A12 ) 



T(n + l/2)i ^ 
g-^ui/^fy /t^woy 

n=0 

Repeating the above sequence P times allows us to synthesize the net transformation 

AT+l 



W aM = (w aM1 w aM2 ) p « J [ e ^«(^K./a i (A13) 



/ 7r(n + l/2)Z > 

n=0 

Where <i# fc ; = a k i/2P. For P sufficiently large, we have ~ W 7 ^. Since this shows we 

can generate arbitrary terms in the Fourier expansion of the Euler angles which define the 
rotation in Eq. (14), we have proven that we can synthesize arbitrary transformations of the 
form, Eq. (14). 

2. Relating P and Q to errors synthesizing 

In this section, we relate the error in synthesizing to P and Q. This relation will 
allow us to discuss the trade off between the accuracy and time with which we synthesize 
in the next section. Direct calculations show us that if 

N+l 

T = Y[ e - in<t>a ^ /2 . (A14) 

n=0 

then the error in synthesizing T a is 

e' T =\\Pc(T-T a )P c \\ (A15a) 
^^'^f (A15b) 

VQ 
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Where we have neglected terms of order ^q 1 ^ and smaller. 
Furthermore, the error in synthesizing any particular W k \ is 



~-w kl 



Pc (W kl - W aM ) P C \\<P\ \Pc (Wl' p - w aM1 w aM2 ) P c \\. (A16) 



If we define 



N+l 



N + Vg m 

J~J e -ide fciv ^(cos( ^ ( ^( 2) ' )^, n -sin( ^ ( '] v h t L ( 2) ' )g- i) , n )/2 

In d9 kl 



n=0 

Wkl2 =TU 2 [ g ma x, 0, 

JV+1 



N + l'g n 



J~J e -^e fc i V ^(cos( ^ ( ^( 2)i )ff J ,, n +sin( ^ ( '^ 4 L ( 2) ' )o- i) , n )/ 



?1=0 

then the error in w a ,kij is 

e 'w aMj = \\Pc{w klj -w aMj )P c \ 
we can bound e' u 



Since e'„ ~ e'„ 

w a,kll w a,kl2 



~-w kl 



^(ll^(< p -w^)p c || + 26; Qfc(i ). 



Furthermore, 



and 



Pc (wl' p - w m w kl2 } P c \ | < V2 {^{N + 1)) 
Where we have neglected terms of order ^ N+ ps ' and smaller. Then 



4-„ < p 



(A17a) 



(Al7b) 

(A18) 

(A19) 

(A20) 
(A21) 



(A22) 



Now we relate the errors in the individual Fourier terms of the Euler angles, W a>k u to the 
errors m Ua , 

£'(2) < 3 (7v + 1) max 1 1 P c (W kl - W kl , a ) P c \\. (A23) 

k,l 

Since the bound on \\Pc (W k i — W k ^ a ) Pc ||, Eq. (A22), is independent of k and I, we have 



E' {2) < 3P V2(N + 1 



2tt 



16 



Q 



(A24) 
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3. Time to synthesize 



Let T(M) be the time to synthesize a transformation M, then 

T(UW) < 3(A + 1) max T(W aM ). (A25) 

kl 

From Eq. (A13) and Eq. (All) we see that 

T(W aM )=4PT(T a ). (A26) 

The time to synthesize T a is bounded by T(T a ) < 4T gX Jl^ so that T(W aM ) = -^Lp^/QT g , 

where T = — . Then 

9 9 

T(tf < 2 >) < -^PVQT g (N + 1). (A27) 
V 2tt 

We are free to choose P and Q, but comparing Eq. (A24) and Eq. (A27), we see there is a 
trade off between minimizing E'^ and T(U^). For most cases, we have a target accuracy, 
E'( 2 \ and we would like to minimize the time required to reach that accuracy. We do so by 
first choosing P and Q so that equality holds in Eq. (A24), which gives 

18432,3 (iv + 1)7 p 4 
(12a/2tt 2 (JV + 1) 3 -£'( 2 )P) 

Plugging this value for Q into Eq. (A27) and finding the value of P which minimizes 
gives 

T(U^) < kVT g {N + V™ . (A29) 

Where w 2.7 x 10 9 . 

There are several reason why the scaling with and E^ is so bad. The first problems 
arise from the approximation used to synthesize t a . Even though t a is effectively a very 
small rotation, the time to synthesize it is still proportional to 1/ y/Q. Then t a is repeated 
Q times to synthesize T a and in order for the approximation Eq. (A9) to hold, Q must be 
large, and so synthesizing T a requires a time much larger than T g . To make matters worse, 
T a is used to synthesize the effectively infinitesimal pulses, Wkij t a- Because Wkij, a is repeated 
P times to synthesize Wki, a , we must repeat T a . In order for the approximation in Eq. (A12) 
to hold, P must be large and because we use T a so many times, it must be synthesized 
with extremely high precision. In addition, because T(wkij, a ) is finite and Wkij, a is repeated 
P times, T(W k ^ a ) is extremely large, which in turn forces our estimate for T(U^) to be 
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extremely large. Thus, the combination of the approximations in Eq. (A9) and Eq. (A12), 
and the finite time required to synthesize t a conspire to generate terrible scaling with N and 
E^ a \ We conclude by noting that vffy is a special case of U^ 2 \ so we have also shown in this 

(2) 

section that we can synthesize u^. 



Appendix B: Synthesis of ufy- 

Next, we construct transformations of the form Eq. (25a). We begin by noting that we 
have already shown we can synthesize transformations of the form Eq. (26a), because they 
have the same form as Eq. (14). Thus to construct ufy we can use either Eq. (27a) or 
Eq. (31a), depending on whether the torque vector of has a z-component. 

the bound on T (v^j-^ is given by Eq. (A29) we have 



The time required to synthesize ufy a is dominated by the time to synthesize V^j- a . Since 



r(.#, J < 2*« 2 >r„^±ip. (Bl) 

However, we would prefer to relate this to the accuracy with which we synthesize vfy- , e'^. 
Since the errors in synthesizing a are due entirely to the errors in synthesizing V^j- a , we 
have e'« = 2E'^, and so 

nuX) £ ^^P- (B2) 

Where fcW « 4.4 x 10 10 . 



Appendix C: Treating the overall phase 

In many cases, we are not concerned about the overall phase of the transformation we 
seek to synthesize and so we are usually interested in minimizing 

e (a) = min \\P C (U (a) - e^U^) P c \ I . (CI) 

4> 

With a = 1,2. However neither the error in T a or W k i is improved by ignoring the overall 
phase, 

min \\P C (T- e*T a ) P C \\=\\P C {T- T a ) P c \\ , (C2a) 
min P c (Wl' p - e^w kll w kl2 ) P c 1 1 = 1 1 P c (w^ /p - w kll w kl2 ) P c \\. (C2b) 
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As a result, ignoring the overall phase does not improve the error or shorten the time to 
synthesize a pulse when using the analytic construction of the appendix. Thus, we can 
replace e'^ with and the calculated times remain unchanged. For comparison with the 
numerical optimization presented later we give the time to synthesize a pulse in terms of 



e («) 



T(UW) <fc ( V%) 3 )10 ' 5 . ( C3a ) 
T{u% a ) <k^T g ^±^. (C3b) 



Appendix D: Time to synthesize arbitrary transformations 

Now we are ready to give the time required to synthesize an arbitrary unitary transfor- 
mation on the computational space. In Sec. Ill we described how to synthesize an arbitrary 
transformation U using ufy. Since we now have a method to synthesize ufflf ~ v>j/!r we can 
use u^jj\f a to synthesize a transformation U a that approximates our target transformation, 
U, with an error bounded by 5, 

mm \\P C (U - e^Ua) P c \\ < 5. (Dl) 

As described in App. B, ufy- a is synthesized using several tx^ Q 's. As a result, all the errors 
which result in synthesizing U a are due to the errors in synthesizing u^ a . If g a (N) = 
3N 2 + 5N + 2 is the number of vfflf used to construct U a , then the error with which we 

(2) 

synthesize any given Ujjj- a must be bounded by 

(2) 



Then the total time for one of the u^ 's used to synthesize U a is 

n^j <- ^ (3Ar2+5Ar+ / (Ar + 1) ' a5 . m 

(2) 

Finally, because the time to synthesize U a is dominated by the time to synthesize and 
there are g a (N) such transformations needed to synthesize U a , the total time is 

nu.) < ^ MLtjM^jg . (D4) 
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It's also common to deal with a quantity which can be related to the gate fidelity, [40], 



''I = m J n mJ, n \\Pc (U - e<*U„) P c \\ 2 . (D5) 



In terms of r\ the total time is 



T(TT \ < hT (3iV 2 + 5iV + 2) 4 (iV + l)9 

T([/ ) < kT g -p . (D6) 



Where jfe « 3.4 x 10 8 . 
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